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Abstract 

A set of simple exactly solvable potentials are shown to have convergent WKB 
series. The resulting all-orders quantisation conditions provide a unified description 
of all known cases where an exact WKB quantisation condition has been obtained by 
modifying the potential (a la Langer), together with several new examples. 
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Investigations, both analytic [1] and numerical [2], into the higher-order behaviour 
of the series (in Ti 2 ) of corrections to the WKB quantisation condition for particu- 
lar potentials have usually shown the series to be divergent. This is in accord with 
the belief, on the basis of general properties of differential equations, that the se- 
ries is generically divergent, though still asymptotic to the correct result [3]. On 
the other hand, there is the long-known fact that the harmonic oscillator and Morse 
potentials have WKB series with exemplary convergence: the lowest order quanti- 
sation condition is exact and all the corrections are identically zero. Additionally, 
Bender, Olaussen and Wang [2] have pointed out that the Rosen-Morse potential 
V(x) = — Vq/ cosh 2 x has a series with a finite radius of convergence, whose terms can 
be calculated and then summed to give a simple all-orders quantisation condition. 
Because the Rosen-Morse potential is an exactly solvable one, it can be checked that 
the summed series gives the correct spectrum, i.e. that no non-perturbative contri- 
butions such as terms like e~ x l n have been missed. The main purpose of this paper is 
to show that this convergent WKB series is not an isolated example amongst simple 
exactly solvable potentials. 

We approach the problem from a slightly different angle than Bender et al by 
noting that the all-orders condition for Rosen-Morse could have been deduced directly 
from the already known spectrum. In general such a quantisation condition will be 
of the form 



where the function f(h, E) contains the o(h 2 ) and higher corrections to the familiar 
WKB condition. Provided one knows E(n) exactly to begin with and can evaluate 
the integral on the left hand side exactly, the only unknown is /. Deriving it this 
way automatically includes any non-perturbative terms. Any possible query about 
the WKB series for V(x) can be answered once / is known. 

Unfortunately, while plenty of exactly solvable potentials are known, V(x) is usu- 
ally sufficiently complicated that the integral is forbidding. However it has proved 
possible to carry out this programme for a set of simple, well-known potentials in- 
cluding the Coulomb, Eckart and Poschl- Teller ones. Introducing a function u(x) 
such that V = u 2 (having fixed V = as the bottom of the well), this set can be 
characterised as containing all the solutions to either 




(1) 




(2) 



or 



— = a + bu 2 + cuv 'a + bu 2 , 
ax 



(3) 
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where a, b and c are arbitary constants. We should point out immediately that there 
is a close connection here to the potentials discussed in [4] in connection with exact 
lowest-order SWKB quantisation conditions (see below). There equations (2) and (3) 
were satisfied by a superpotential <j>, giving potentials 

V(x) = 4> 2 - h<p' + e = u 2 , (4) 

where eo is the ground-state energy. As detailed in the Appendix, any solution u(x) 
to either (2) or (3) with constants a, b and c corresponds to a potential V(x) with a 
(f>(x) satisfying the same equation with some constants A, B and C . Thus the set of 
potentials considered here is exactly the same set discussed in [4]; while not essential 
to the argument, the fact that these potentials have exact SWKB conditions greatly 
simplifies the algebra. A minor exception aside [5], these potentials are usefully 
tabulated by Dutt, Khare and Sukhatme [6].[] 

Note that equation (2) includes as solutions both the trivial harmonic oscillator 
and Morse examples (for which / = 0) and the Rosen-Morse example of Bender et 
al. This paper will thus directly generalise these earlier results. 

We begin by proving perturbatively that for this set of potentials / is independent 
of E. Of course this will become quite plain when we derive / non-perturbatively, 
but the proof does bring out some of the significance of this set for the SWKB 
approximation. Making the substitution if) = exp(iS/h), the Schrodinger equation 
can be written as (2m = 1) 

S' 2 - ihS" + V = E (5) 

and the WKB approach involves expanding S as a power series in fi, the coefficients 
of which can be generated recursively [8]. Once connection formulae have been con- 
sidered, one obtains a quantisation condition which is a series in h 2 ; in terms of 
u(x) 

The higher-order corrections quickly increase in complexity, but any term appearing 
in the o(h 2n ) one must have the structure 

d M f u^u^ 
dE™ J (E — m 2 ) 1 / 2 dx 

1 These potentials all share the elegant property of "shape invariance" (reviewed in [6]) and the 
emphasis in [4] was on the relation between it and the SWKB series. Recent developments [7] have 
weakened this relation, but do not effect the conclusions in [4] concerning potentials for which the 
lowest-order SWKB condition is exact. 
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ni + n 2 + . . . + n m = 2n, M = (2n + m - 2)/2 

u(°)=ii, u ( 2 ) = «",.-■ (7) 

It can also be arranged that n m = 1, so the integration variable can be changed from 
x to u. Now notice that equations (2) and (3) have the unique property that 

d N v 

^ = p(N + l)+p(N)g(u), (8) 

where p(m) is used to denote an otherwise unspecified rath order purely odd or even 
polynomial in u. Solutions to (2) have g(u) = 1 and those to (3) g{u) = \/a + bu 2 , 
but no other choices of g repeat this pattern (8). For these cases, (7) becomes 

d M [Ve p (2M) +p{2M- l)g{u) 



dE M J-^E (E - u 2 ) 1 / 2 

The term containing g(u) is odd and integrates to zero; the other is a polynomial in 
u 2 with leading term u 2M , so that the integral gives a polynomial in E with leading 
term E M . Differentiating M times leaves a constant independent of E. Since this 
happens in all terms in the series, for these potentials f(%, E) reduces to f{%) with 
no -E-dependence. Their all-orders quantisation condition (1) is thus particularly 
similar to the standard lowest-order one: the .E-dependence is exactly the same and 
the higher-orders can be treated as an additional constant on the right hand side, 
absorbable into a shifted ra[| There may be other potentials with ^-independent 
higher-orders, but experience with the SWKB series suggests otherwise [4]. 

Note that the SWKB series would be derived by splitting V into <p 2 — h<p' and 
treating the second piece as o(H) [10]. This leads to a slightly different approximation 
scheme with a series somewhat like (6), but with integrals over rather than u. The 
arguments above go through rather similarly, except that the presence of a factor of 
E — eo multipying each correction (the SWKB approximation is exact as E — > e ) 
means that the final integral is over a polynomial p(2M — 2) and thus that the M 
derivatives leave nothing. Each correction is zero and these potentials have exact 
lowest-order SWKB conditions [4] 

J yjE-e - <f) 2 dx = rmh. (10) 



2 Potentials with / independent of E are the only ones for which the semi-classical inverse method 
[9] is exact. This is because differentiating (1) with respect to E (the first step therein) then gives 
a formula indistinguishable from the lowest-order one. Given n(E), one can thus reconstruct the 
"excursion" of V(x) exactly. 
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It is this fact that allows for a particularly simple derivation of / for these poten- 
tials, avoiding the necessity of calculating it separately for each individual potential. 
Because of (10), we can rewrite (1) as 



/ 



VE - u 2 f VE-e o -0 2 irk 

^^du = j dcf> + - + nf. (11) 



But we know u' in terms of u from (2) and (3) and similarly for <//, so the integrals 
are seen to be elementary ones. That on the right gives a function of A, B and C, 
but this can be recast in terms of a, b and c using the formulae in the Appendix. For 
potentials corresponding to (2) the result is 



As expected, the ^-dependence has cancelled and / — > as b — > 0, the latter to 
give agreement with the known absence of higher-order corrections for the harmonic 
oscillator and Morse potentials. Setting b = l/\Ao shows (12) to be in agreement 
with the Bender et al result [2]. 

The equivalent function for potentials given by (3) is more complicated. 

H 2 (h 2 ) = I (46" 1 (6+c 2 ) + h 2 (c 2 - bf - (c 2 - 6)^/(4/6 + h 2 b + h 2 c 2 ) 2 - Ah 4 c 2 b) . (13) 
8 ^ ' 

Note that H(h 2 ) — > 1 as c 2 — > 6, so that, contrary to appearances, the c 2 = 6 

case is not pathological; this corresponds to the radial equation for the 3D harmonic 

oscillator and has 

/ = l(vT7ftV-i). (14) 

These last equations, (12)-(14), are the central results of this paper. All these 
functions are analytic about H — 0. 

Equation (12) is particularly important in that it sheds new light on a classic body 
of work concerned with modifying potentials such that the lowest order WKB condi- 
tion for this new potential gives the correct spectrum for the unmodified one. This 
can clearly be done in principle for any potential, but a surprisingly simple modifi- 
cation was discovered for the radial version of the Coulomb potential very early [11]: 
replace 1(1 + 1) in the centrifugal barrier by (/ + 1/2) 2 . Langer [12] subsequently 
showed that a consistent treatment of the WKB approximation for radial problems 
introduces an additional term in the lowest order condition that is just equivalent 
to this modification. But this is a demonstration that the condition for which the 
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Coulomb potential is exact is the correct lowest order one rather than an explanation 
for this exactness. The first proof (excluding empirical observation that the trick 
works) was by Rosenzweig and Krieger [13] within the Froman version [14] of the ap- 
proximation. However this was only as one of five case-by-case examples of potentials 
which could be modified to give the correct answer, some of which were unrelated to 
Langer's work, e.g. the Rosen-Morse potential V(x) = —V / cosh 2 a; can be modified 
using [15] 

V'o = V Q + ^, (15) 

but is an innocuous ID problem. 

Surprisingly, all the potentials considered by Rosenzweig and Krieger are members 
of the set considered above and so it is possible to bring together their stray examples 
into one generalisation. Before turning to this, we should however deal with the 
possible objection that Langer's observations invalidate the derivations above for 
particular cases, notably the Coulomb one. For instance, the all-orders result appears 
to imply a lowest-order quantisation condition in disagreement with Langer's for the 
Coulomb problem. Actually, because Langer's modified potential is invariably H- 
dependent, there need not necessarily be any contradiction between the two lowest- 
order formulae; the modified condition may result from a partial summation of higher- 
order pieces in the original series (6). And for the set of potentials here that series 
(6) can always be rendered meaningful by appealing to their exact SWKB conditions, 
since the WKB series (6) can always be obtained by rearranging the SWKB one (c.f. 
[16]). Thus naively applying (6) to the Coulomb potential gives a series in h which 
will sum to the correct answer and uniqueness requires that this be the same series 
implied in (12). Of course no such justification is necessary for potentials like the 
Rosen-Morse one which have no singularities. 

Rewrite (2) as 

u' — b(u + p)(u + q). (16) 

Hitherto the arbitrary reference energy has been chosen such that V = at the 
minimum of the potential, but this minimum is one of the things that can be expected 
to be modified, so we introduce E = E — p 2 . For the Coulomb case this corresponds 
to refering all energies to the top of the well. From (1) and (12), the exact condition 
is 



-1 

b ' b(p-q) 



- q \[l - (E + p 2 )/q 2 - py /l-(E + p 2 )/p 2 ] = (n + l/2)nh + J 1 + 



1 ?! 

b 2 + 4 
17) 
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and we seek modified b' , p' and q' such that 
-1 1 

T7 + 



V b'(p' - q') 
Thus 



q'Jl-(E + p' 2 )/q' 2 - P yi-(E + p' 2 )/p' 2 = (n + l/2)nh. (18) 



1 1 Ti 2 

772 = 77 + - 



b' 2 b 2 4 

q 2 -p 2 = q' 2 -p' 2 , b{p-q) = b'{p'-q'). (19) 



In the Coulomb case 



-e 2 



b = p = q = == (20) 

Hyjl(l + 1) 2tiyjl(l + 1) 

and, as expected, (19) implies that the substitution 1(1 + 1) — > (Z + 1/2) 2 gives 6', p' 
and </. As a example of a new result obtained, consider 

"^0 i o /tv~x+„„u^ i x2 



V(x) = f- + 2JVb5 tanh x + <5 2 , (21) 

cosh a; v 

a more general version of the Rosen-Morse potential dealt with by Rosenzweig and 
Krieger [13] and Bender et al [2]. The modified potential is given by (15), now 
supplemented with 

* = J¥f (22) 

All the modifications tabulated in [14] are subsumed into (19), along with several 
such generalisations. 

To summarise, we have shown that potentials V = u 2 given by solutions to (2) and 
(3) have - in addition to exact lowest-order SWKB conditions (10) - WKB series (6) 
which have finite radii of convergence and which sum to give the particularly simple 
quantisation conditions corresponding to (12) and (13) respectively. For the potentials 
with singularities, Langer's modification of the lowest-order WKB formula has been 
seen to be equivalent to the summation of the higher-order corrections in the original 
(and here ostensibly incorrect) version of the approximation. In addition, equations 
(19) provide a recipe for modifying the potentials - whether containing singularities 
or not - in all those cases where such simple modifications have given the correct 
spectrum in lowest-order WKB, together with some new examples. 

This sample of potentials is obviously very special and the immediate challenge is 
clearly to find f(h, E) for other exactly solvable potentials. 



Appendix 
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Consider two functions u(x) and <j>(x) related by 

<f) 2 - fifi + e = u 2 . (Al) 

If is a solution to 

4>' = A + B<f) 2 + Ccf), (A2) 

then u is a solution to (2) with 

„-M .mi/2,4 | ^ 2 ( 2 ~^) 

b = Bil-hB)- 1 / 2 , c=C{l-hB)- 1 

Similarly, if is a solution to 

4>' = A + B<f) 2 + + 50 2 , (A4) 

then u is a solution to (3) with 

a = ay/AB - (3CVA, b = a(a 2 - p 2 )' 1 
c = b- 1/2 (pVAB-aCVA) 

eo = hA+ A{l ~ B hB) (l " y/l-h 2 C*B/(l-hB)*), (A5) 

where a 2 = A/ B — e and j3 2 = e — ftA 

In both cases the appropriate converse is also true. 
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